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The dynamics in the thin boundary layers of temperature and velocity is the key to a deeper 
understanding of turbulent transport of heat and momentum in thermal convection. The velocity 
gradient at the hot and cold plates of a Rayleigh-Benard convection cell forms the two-dimensional 
skin friction field and is related to the formation of thermal plumes in the respective boundary 
layers. Our analysis is based on a direct numerical simulation of Rayleigh-Benard convection in a 
closed cylindrical cell of aspect ratio T = 1 and focused on the critical points of the skin friction 
field. We identify triplets of critical points, which are composed of two unstable nodes and a saddle 
between them, as the characteristic building block of the skin friction field. Isolated triplets as well 
as networks of triplets are detected. The majority of the ridges of line-like thermal plumes coincide 
with the unstable manifolds of the saddles. From a dynamical Lagrangian perspective, thermal 
plumes are formed together with an attractive hyperbolic Lagrangian Coherent Structure of the 
skin friction field. We also discuss the differences from the skin friction field in turbulent channel 
flows from the perspective of the Poincare-Hopf index theorem for two-dimensional vector fields. 

PACS numbers: 47.27.N-, 47.55.pb 


I. INTRODUCTION 

One way to gain a deeper understanding of turbulence is to break the complex fluid motion down into a few basic 
building blocks which are connected and often appear repeatedly on different scales. Their geometry, dynamics and 
interaction determines the global statistical properties of the turbulent flow mm- They can be identified either in the 
Lagrangian or Eulerian picture. For example, the most persistent attracting and repelling Lagrangian manifolds in a 
turbulent flow, which are known as Lagrangian Coherent Structures (LCS) [3], are found to serve as barriers to the 
mixing of scalar quantities and thus affect the global scalar fluctuations HU. The stretching of simple Burgers-type 
vortices is used to describe the small-scale intermittency and the related evolution of high-amplitude enstrophy events 
as well as to explain the non-local energy transfer in the turbulent cascade [zhtq!. When a passive scalar increment is 
conditioned to the distance of zero-scalar-gradient points, scaling properties similar to scalar structure functions follow 
mnn]. In all the above examples, geometric objects are connected to dynamical processes and statistical properties 
of the corresponding turbulent flow. 

In this work, we apply these ideas to the formation of thermal plumes in turbulent Rayleigh-Benard convection 
(RBC) which is present in a fluid layer between two horizontal parallel plates heated from below and cooled from above 
m- Thermal plumes are fragments of the thin thermal boundary layer. They are formed in the vicinity of the heated 
(cooled) plate, then rise (fall) into the bulk of the turbulent convection layer and provide the major contribution to the 
near-wall turbulent heat transfer and its local fluctuations. Their dynamics and size have been studied, for example, in 
experiments [HHni and direct numerical simulations [MS]. Here, we will describe the formation of thermal plumes 
by the vertical temperature derivatives in connection to the two-dimensional velocity gradient vector field directly at 
the heated (cooled) plate. It is considered as a blueprint of the coupled temperature-velocity dynamics above the heated 
(cooled) plate. The structure of the velocity gradient field at the plate is determined by its critical or zero points and 
their local topology. As we will see, triplets of such critical points play a central role in the near-wall dynamics and can 
be connected to the formation of line-like thermal plumes. Their dynamics will be analyzed in the Eulerian as well as in 
the Lagrangian frames of reference. We demonstrate that ridges of line-like plumes coincide with attracting hyperbolic 
LCS of the skin friction field. 

The stem of thermal plumes at the hot bottom plate coincides with local minima of the absolute vertical temperature 
derivative \dT/dz\ at the plates, where z is the vertical direction. We will therefore derive and discuss a reduced 
transport equation for this vertical temperature derivative and its magnitude which is based on the fact that at the 
wall, the velocity gradient tensor = duijdxj has only two non-vanishing components which form a two-dimensional 
vector field. We then relate the resulting network of critical points - saddles and nodes [52] - to regions on the plate 
where the thermal plumes are lumped together before their take-off into the bulk. In other words, the eigenvalues which 
determine the phase portrait of the critical points appear in source-sink terms of this reduced transport equation for 
dTjdz. 

We use the notion which has been suggested in Ref. [53] and term the two-dimensional vector field, which is composed 
of the wall-normal derivatives of the tangential velocity components, as the skin friction field. Critical points of the 
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FIG. 1: (Color online) Critical points of a two-dimensional vector field. Phase portraits are displayed, (a) Saddle, (b) Unstable 
node, (c) Stable node, (d) Unstable focus, (e) Stable focus. 


skin friction field have been studied recently in plane shear flows with an emphasis to understand their connection 
to near-wall backflows and their role in the formation of near-wall coherent streamwise vortices, both in simulations 
|23II25] and laboratory experiments [261 El]. It has been reported that critical points are partly found at the tail of 
an evolving large-scale structure [23]. As we will see in the following, there are differences from the turbulent shear 
flow case that has a well-defined unidirectional mean flow. On one hand, the differences arise from the different flow 
geometry. On the other hand, they also arise from the active coupling of the temperature and velocity fields which is 
absent in a turbulent channel flow. 

Our present work aims at applying this framework to the convection flow and plume formation. Here, we outline 
therefore the concepts rather than provide an extensive parametric study with respect to Rayleigh and Prandtl number, 
two dimensionless parameters that characterize convective turbulence. Such a study is certainly necessary and will be 
done elsewhere. This means that the concepts are tested with one turbulent convection data set at a moderate Rayleigh 
number. The analysis in the present work is done mostly at the bottom plate of the convection cell. The outline of the 
paper is as follows. We introduce in brief the Boussinesq model of convection. Afterwards, we lay out the theoretical 
basis and develop a reduced equation for the evolution of the vertical temperature derivative. Our predictions are tested 
afterwards using high-resolution data sets at a moderate Rayleigh number. We close the paper with conclusions and a 
brief outlook. 


II. EQUATIONS OF MOTION 

We solve the three-dimensional Boussinesq equations for turbulent RBC in a cylindrical cell of height H and diameter 
d. The equations for the velocity field u{x,t) and the temperature field T{x,t) are given by 

V-M = 0, (1) 

d'UL 1 

+ {u- V)u = -Vp -b vV'^u + ga(T - To)e^ , (2) 

at po 

rfr 

— ^{xi-V)T = kV^T. (3) 

at 

The pressure field is denoted by p{x,t), the constant mass density by po and the reference temperature by Tq pOl 
The aspect ratio of the convection cell is T = d/H = 1. The Prandtl number which relates kinematic viscosity ly and 
thermal diffusivity k is Pr = vj k = 0.7. The Rayleigh number is Ra = gaATH^/{i/K) = 10^. The variables g and a 
denote the acceleration due to gravity and the thermal expansion coefficient, respectively. The temperature difference 
between the bottom and top plates is AT. In a dimensionless form all length scales are expressed in units of i7, all 
velocities in units of the free-fall velocity Uf = yjgaATH and all temperatures in units of AT. We will denote the 
heated plate by M. 

We apply a spectral element method in the simulations in order to resolve the gradients of velocity and temperature 
accurately |30j . More details on the numerical scheme and the appropriate grid resolutions can be found in Ref. uniiis]. 
No-slip boundary conditions are applied for the velocity at all the walls. The top and bottom walls are isothermal and 
the side wall is thermally insulated. The cylindrical convection cell is covered by 30720 spectral elements. On each 
element all turbulent fields are expanded by llth-order Lagrangian interpolation polynomials with respect to each 
spatial direction. 













3 


III. THEORETICAL BASIS 


A. Skin friction and near-wall velocity fields 


At the bottom plate (2 = 0), the velocity gradient tensor Aij takes the following form 

/O 0 
A 


2=0 


0 0 duxjdz 

= I 0 0 duyjdz 
0 0 0 


Both components form a two-dimensional wall shear stress vector field which is defined as 

9m j 


'^w — 


dz 


(4) 


(5) 


2 = 0 


where the subscript _L denotes the two horizontal (or tangential) x- and y-components or coordinates. Skin friction 
and wall stress fields are dehned by derivatives normal to the plates. Therefore df jdz has to be changed to —dfjdz 
for any held / when the top plate behavior is analyzed. The skin friction held is dehned by 


s — 


Pov 

Note that the vorticity vector at the wall, which is denoted as surface vorticity [53], is given by 

duy dua 

wt = I - 


( 6 ) 


dz' dz ' ’ 

and is perpendicular to s. In the vicinity of the top and bottom walls, we expand the velocity components as follows: 

Mj.(a;_L,2,t) = s(a;_L,t)z, (8) 

(9) 


1 


M2(a:_L,2,t) = • s)z‘^ + C'(a;_L). 

Again, = {x,y). The integration constant C{x±) = 0 in order to satisfy the no-slip boundary conditions at the 
wall. The velocity held which is given by Eqns. ([8) and ([^ satishes the incompressibility condition. It can be already 
seen that the essential information of the near-wall velocity is captured by the skin friction held, in particular by its 
divergence. Sources and sinks of s are connected to the critical points. 


B. Critical points of the skin friction field 


1. Classification 


The local topology or phase portrait of the critical points of the skin friction held is determined by the eigenvalues 
of the corresponding Jacobian which is given by 


J = 


characteristic equation is solved: 

• Saddle: Ai = ai < 0 and A 2 = 02 > 0 

• Unstable Node: Ai = oi > 0 and A 2 = 02 > 0 

• Stable Node: Ai = oi < 0 and A 2 = 02 < 0 


/ d'^Ux 

d'^Ux \ 

dxdz 

dydz 

d'^Uy 

d'^Uy 

\ dxdz 

dydz / 

pairs of 

complex 


( 10 ) 


• Unstable Focus: Ai ^2 = a±ib (a > 0) 

• Stable Focus: Ai ^2 = —a ±ib (a > 0) 

The local topology of all these critical points is sketched in Fig. Stable nodes or foci can be considered as sinks while 
unstable nodes and foci are sources. 
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2. Poincare-Hopf theorem 

The Poincare-Hopf theorem states that for a vector field with isolated zeros that is defined on a compact orientable 
differentiable manifold the Euler characteristic x(A4) can be related to the sum of the indices of the N isolated 
critical points x* [ST], 


N 


^ind(a;*) = x{M). 


k=l 


Here, the index of a critical point x* in a two-dimensional plane is defined as, 


ind(a:*) = 


1 

27r 


dO , with 9 = arctan 


C(x*) 


( 11 ) 


( 12 ) 


Curve C{x*) surrounds the critical point. For a stable or unstable node and a focus, we find ind(a;*) = -|-1, and for a 
saddle ind(a;*) = — 1. The Euler characteristic is a number that can be assigned to the topological properties of the 
plane, which is the two-dimensional manifold in the following. In the wake of this theorem, a first significant difference 
between turbulent channel flows and turbulent convection cells arises for simple geometric reasons. DNS of a turbulent 
channel flow use typically periodic boundary conditions in both horizontal directions. The no-slip walls of the channel 
are thus two-dimensional tori and have an Euler characteristic of zero. Consequently, Eq. (11) implies that the number 
of nodes (and foci) of the skin friction field equals the number of saddles [231 [32] ■ This is different to the present 
convection case. The bottom and top plates M are orientable circular disks with a boundary and hence the Euler 
characteristic is x(M) = 1. We will come back to this point when we analyze our simulation data. 


C. Hyperbolic Lagrangian Coherent Structures 

As mentioned before, we denote the bottom plate of the convection cell by M. The time-dependent skin friction 
field s generates a flow map : M ^ M which advects Lagrangian trajectories from time to to t = to + r, i.e., 

^to^^(®o) : I—>■ x{t',Xo,to) in correspondence with 


X = s{x, t). 


(13) 


The Cauchy-Green strain tensor is given by 


C(xo) 




(xo) 


^^0+r 


(a;o), 


(14) 


where the asterisk denotes the transposed matrix. The finite-time Lyapunov exponent (FTLE) follows as 

A^(a;o) = ^ log CTr with (Tr = Xniax{C{xo)) . (15) 

ZT 

Lagrangian coherent structures are the most attracting or repelling hyperbolic manifolds over a finite time interval [3]. 
In typical situations, ridges in the FTLE field indicate the location of LCS [331, so our approximation will be based on 
this. In contrast to many previous applications of LCS to fluid dynamical problems, here we will examine the structures 
in the two-dimensional vector field s that is not divergence free. The transport of temperature in the velocity field in 
the vicinity of the bottom plate, i.e., at z 1, is translated to a transport of the vertical temperature derivative in the 
skin friction field at z = 0. The corresponding transport equation is derived and simplified in the following section. 


IV. TRANSPORT OF WALL-NORMAL TEMPERATURE DERIVATIVE 


We denote the temperature gradient by G = {G±, Gz)- In the following we will derive an equation of motion for the 
vertical temperature derivative, Gz = dT/dz, very close to the bottom plate at z = 0. This quantity can be directly 
related to the plumes (see also Fig. [plater in the text). Line-like thermal plumes are formed where \Gz\ is very small, 
whereas regions with enhanced velocity shear will exist for large \Gz\, i.e., \Gz \ S> 1- The equation for Gz follows from 
Eq. (§. Evolution equations for scalar derivatives have been discussed for example in Refs. [331IM] in the context of 
passive scalar mixing. The evolution equation for the vertical temperature derivative reads 


dGz 


+ {u-V)Gz 


du 

dz 


G + kV^Gz- 


dt 


(16) 
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With the near-wall expansions the transport equation for the vertical temperature derivative ( |16| ) close to the wall 
can be rewritten as 


dG, 

dt 


s ■ [zV^G, + Gj.] = (Vi • s) 


2 dz 


G. 


2 

Vi H- 

^ dz^ 


G.. 


(17) 


We make several approximations near the wall and thus leave the rigorous treatment. In the spirit of an asymptotic 
expansion, we compare the magnitude of the different terms in (17). Therefore, we set x = xH, y = yH and z = ezH 
where the tilde stands for dimensionless quantities and e <C 1. Thus dx and dy ~ 1/H and dz ~ \/{eH). As a 
consequence, the horizontal contribution to the Laplacian can be neglected in 0-' Furthermore, the two contributions 
to the first term on the right hand side are of the same order of magnitude and can be summarized to one term. Close 
to the wall, we can additionally assume that [G^l |Gj^|. Together with the e expansions, this simplifies (17) to 


dGz 

dt 


sz ■ V^Gz « (Vj^ • s)zGz 


d^Gz 


(18) 


Multiplying the equation with Gz results in the following balance equation for the vertical temperature gradient mag¬ 
nitude 


dGl 

dt 


2(V_l • s)zGI + 2kGz 


d^Gz 

dz^ 


(19) 


The total time derivative on the left hand side of (191 summarizes the partial time derivative and the advection which is 
dominantly by the horizontal velocity components. If we align the coordinate axes in the vicinity of the critical points 
with the two eigenvectors, i.e., {x = 1) and {y = 2), Eq. (19) can be further simplified to 


dGl 

dt 


2(Ai -l- A 2 ) zG^ 2 kGz 


d^Gz 

dz^ 


( 20 ) 


The first term on the right hand side of ( |20[ ) is a gain term when unstable nodes or foci are locally present. In turn, 
the magnitude of Gz decreases when stable nodes or foci exist. This simplified balance equation illustrates clearly the 
close connection between the local topology of the skin friction field which is given by the eigenvalues Ai and A 2 and 
the vertical temperature derivative and thus the plumes. 

As a short note, we briefly discuss the case of free-slip boundary conditions at the heating and cooling plates which 
are given by 


Uz{x,y,z = 0) = 0 


and 


dux 


dz 


2=0 


duy 

dz 


= 0 . 

2=0 


( 21 ) 


In this case the skin friction field vanishes, s = 0, and yet the plumes are generated. In such a scenario, the velocity in 
the vicinity of the wall can be expressed to the first order as 


u±{x±,z,t) = Uw{x±,t), ( 22 ) 

Uz{x±,z,t) =-{'V± ■ Ui,)z, (23) 

which are the free-slip counterparts for eqns. ([^ and respectively. Here, = (ux,Uy) represents the velocity 
field on the wall with sources and sinks (see also Ref. [35]). Using the procedure similar to that used for the no-slip 
boundaries, we obtain the approximate evolution equation for Gz which is given by 


dGz 

dt 


(V _L • ri^) Gz ^ 


d^Gz 
dz^ ■ 


(24) 


It is very clear from the above equation that in the case of free-slip conditions, the divergence of the slip velocity 
field Uxu acts as the source for the evolution of Gz as compared to V_l • s (see equation ( [I^ ) in the no-slip case. In 
other words, with free slip boundaries, the two-dimensional slip velocity on the wall plays the role of skin friction field 
in characterizing the plume dynamics. In this paper, we only present and analyse data considering no-slip boundary 
conditions. For a quantitative analysis we will now turn to the direct numerical simulation data. 


V. SIMULATION RESULTS 
A. Skin friction field at top and bottom plates 

Fig.i illustrates the relation between the skin friction field, in particular its critical points, and the vertical tem¬ 
perature derivative for one time snapshot of the simulation. Data for the plane z — 0 are shown in the upper row. 
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FIG. 2: (Color online) Simulation snapshot which displays the critical points of the skin friction field and thermal plumes. Top 
row: z = 0, bottom row: z = 1. Left: Magnitude (colored background) and streamlines of the skin friction field s. Unstable 
nodes and saddle point are indicated by arrows. Middle: Corresponding contour plot of the vertical temperature derivative, 
dT/dz. Line-like plumes are indicated again by arrows. Right: Contours of the divergence of the skin friction field together with 
the streamlines of s. Data are for Ra — 10^, Pr = 0.7 and F = 1. 


data for z = 1 in the lower one. The left contour plots display the critical points of s as the darkest regions of the 
contour background which stands for |s|. The overlaid streamlines of s indicate the local topology at the critical points. 
We highlight a saddle and an unstable node, the latter of which are sources of the skin friction field, at both plates. 
Unstable nodes are generated by cold downward flows which hit the bottom plate or hot upward flows which hit the 
top plate. They are surrounded by regions of high shear as indicated by the magnitude of s. The held lines of the skin 
friction held point almost perfectly radially away from the node, a source region of the skin friction held. 

In the mid panels we display the contours of dT/dz at z = 0 (top row) and —dTjdz at z = 1 (bottom row). Contours 
of dT jdz equal contours of the temperature T slightly above the plate. This was verihed, but is not displayed here. The 
regions with the largest vertical temperature derivative, i.e., with the smallest magnitude, are considered as the regions 
in which line-like plume ridges form and detach subsequently. They coincide with the regions where the magnitude of 
the skin friction held - and thus the tangential stresses or the local shear - are small. High-shear regions fall in between 
regions in which cold plumes hit the bottom plate or hot plumes are about to detach. 

The comparison of the skin friction helds at the bottom and top plates shows also that both helds are not synchronized 
by hot line-like plumes rising from the bottom straight to the top and cold line-like plumes falling from the top down 
to the bottom of the cell. Plumes will be dispersed by the turbulence in the bulk of the cell when they rise and fall. 
The impact of colder fluid at the bottom and hotter fluid at the top, i.e., the process which generates unstable nodes of 
the skin friction held, is thus most probably a consequence of the incompressibility. The rise of plume ridges into the 
bulk pushes neighboring huid back towards the bottom plate. It can be expected that this behavior depends also on 
the Prandtl number of the working huid. 

The lower right plot of Fig. shows the divergence of the skin friction held. Local minima of |s| can coincide with 
unstable nodes. This is however not a necessary condition as the hgure shows. The data support our simple dynamical 
model which is condensed in Eq. ( [20| . Critical points in the form of unstable nodes, for which Ai -|- A 2 > 0, are sources 
for G^, i.e., the magnitude of G^ has a local maximum in their vicinity. In these regions, cold plumes hit the bottom 
plate from above. Sinks of the skin friction held coincide with plume formation regions. When Ai -I- A 2 < 0, the hrst 
term in (201 turns into a loss term for G^. 
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FIG. 3: (Color online) Probability density function of the vertical temperature derivative dTjdz which is taken at the bottom 
plate, z = 0. Data are for Ra — 10^, Pr = 0.7 and F = 1. The dashed vertical line indicates the mean vertical temperature at 
the bottom plate which equals the Nusselt number of the convection flow. 


B. Derivative statistics 

The vertical derivatives of the velocity and temperature are discussed in Fig. In this figure, we display the 
probability density function (PDF) of the vertical temperature derivative at z = 0. It has negative amplitudes only 
and develops an extended tail. Data have been obtained from 230 statistically independent samples of the turbulent 
convection flow which have been gathered along a long-time integration. They are separated by at least one free-fall 
time Tf = H/Uf. Figure |4(a)| display the statistics of the two components of the skin friction field obtained from the 
same dataset. The whole bottom plate is used for the analysis in the following. We verified that the exclusion of side 
wall boundary data does not change the results. The distributions display the typical exponential tails which reflect 
the intermittent statistics of velocity derivatives. 

It is interesting to note that the mean values of the distribution differ (s^ = —3.12 and Sy = —0.49) indicating 
a large scale circulation (LSC) which breaks azimuthal symmetry. This is readily confirmed from the streamlines of 
the time-averaged skin friction field shown for z = 0 in Fig. |4(b)[ that are mostly aligned in the negative x-direction. 
Accordingly, we see an impact region on the right and a “take off“ region on the left. In addition, the left tail of the 
PDF of Sx is seen to be systematically fatter than that of Sy in the range when |sa;| <40 and the opposite ocuring for 
the right tails. It is known from other DNS studies [MISS] that LSC exhibits a very slow angular drift (at time scales 
larger than 1000 free-fall units Tf), which explains the biased directionality of the LSC in our case due to the fact that 
the time window for averaging was smaller. In addition, the LSC depends on Ra and becomes less coherent at large 
Ra. The Rayleigh number of 10^ in our case being relatively small, we thus see a clear indication of the LSC. 

As a next step, we connect the statistics of the vertical temperature derivative with that of the skin friction field. 
Figure shows the statistics of V • s. While one graph displays the unconditioned divergence field, the other graph is 
conditioned to the highest magnitude temperature derivatives, i.e. to events with dT/dz < {dT/dz)A,t where {■)A,t is 
a plane-time average taken at z = 0. The whole positive tail of the conditioned PDF exceeds that of the unconditioned 
PDF. This finding statistically underlines the observation, which was discussed earlier: positive divergence of the skin 
friction field is connected with unstable nodes which hit the bottom plate and generate high magnitude temperature 
derivatives. 


C. Eigenvalue analysis of critical points 

The critical points of the two-dimensional skin friction field are usually not located on the grid points where all 
the discretized field variables are approximated. Due to this reason, a simple search for them ‘on’ the grid points is 
insufficient. To determine the critical points, the following procedure is adopted which is similar to the observation 
that was discussed in Cardesa et al. [32]. At each grid point k, considering the local gradients of s, the position 
dr^ = [dx^, dy*’] (relative to the grid point being considered) at which s = 0 is computed by solving 


j>^dr^ = 


( 25 ) 
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FIG. 4: (Color online) a) Probability density function of the two components of the skin friction field, du^jOz and duyjdz at 
z = 0. The data set for this analysis is the same as in Fig. Mean values are Sx = —3.12 and Sy = —0.49 respectively, b) Lines 
of mean skin friction field averaged over 202 time snapshots which are separated by roughly one free-fall time unit T/ = H/Uf 
of each other. The background contour stands for dTjdz. 



ir _^_ 1_, a, 

-2000 0 2000 4000 

V • s 


FIG. 5: (Color online) Probability density function of the divergence of the skin friction field. We compare the unconditioned 
divergence with the divergence for which dTjdz < {dTjdz)A,t- The data set for this analysis is the same as that in Fig.|^ 


where the superscript k on J and s indicates values computed at the grid point k. Any point so obtained, which is 
farther than 5 from the grid point k is discarded (i.e. if \dr^\ > 6, S being the grid size) and we denote the remaining 
points (which are ‘potential’ critical points) as the set P. 

In the next step, all the points in P are traversed to identify ‘clusters’ of points, where a cluster refers to a group 
of points that are within a distance from each other. In this process, points that do not have at least a single 
neighbor are discarded as outliers. Typically each cluster contains 2 to 4 points which are located within a grid cell and 
were correspondingly obtained from the corner grid points of the cell. Finally, from each of these clusters one point is 
chosen randomly to qualify as a critical point and the others discarded. In this way, it is ensured that all the critical 
points are identified accurately. The type of the critical point is determined from the characteristics of the eigenvalues 
of J as described in subsection III. C. For this analysis, the fields on the unstructured spectral element grid have been 
interpolated to a uniform Cartesian mesh, which simplifies the search of critical points significantly. 

In Fig. (top) we display the results of the eigenvalue analysis using 202 statistically independent time snapshots of 
the DNS. The number of foci is small compared to the number of saddles and nodes. It can be observed that the number 
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FIG. 6: (Color online) Top: Number of saddles, nodes and foci for the individual snapshots of the turbulent convection flow. The 
dashed lines with same color as the data correspond to the arithmetic means. Bottom: Sum of the indices of saddles and nodes. 
The dashed line at 1 indicates the Euler characteristic for the hot bottom or cold top plates M. Same data are taken as in top 
figure. 


of critical points is strongly fluctuating. Furthermore, we detect that the number of nodes and saddles is not the same 
(see bottom panel of Fig.[^. This is a significant difference from the numerical studies in the channel flow. The reasons 
are as follows. First, the Euler characteristic is not zero for our circular disk with boundary as we have discussed in Sec. 
Ill B. Secondly, due to the no-slip boundary conditions, the skin-friction field vanishes on the boundary, i.e., s = 0, 
so that the Poincare-Hopf theorem for manifolds with boundaries is not applicable. Thirdly, for the same reason even 
the critical points in the interior of the disk are often found very close to sidewalls. We note that an extension of the 
Poincare-Hopf theorem m holds also for the case of isolated critical points on the boundary of M. Thus the search 
algorithm for the critical points excludes a very tiny region at the boundary of the circular bottom and top plate, in 
order to restrict to saddles and nodes away from the sidewall. Only these critical points are in the focus of this analysis 
and remain an important structural element of the skin friction field of thermal convection compared to channel flows. 

Figure [T] underlines that the plume ridges can be associated in most instances with node-saddle-node triplets. These 
triplets can be found isolated as in Fig. or they can form a whole triplet network as shown in this figure. Fluid is 
expelled radially outward from the unstable nodes which are distributed over the plate. Centered between two such 
nodes is a saddle point of s which is formed in a stagnation point flow configuration. A similar behavior is observed in 
other snapshots which we have analyzed in detail. It is needless to say that exactly the same dynamics is present at 
the top plate for the formation of cold plumes that fall into the bulk of the cell. 

How can the regions around unstable nodes be separated from each other? Figure [^displays field lines of the surface 
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FIG. 7: (Color online) Network of critical points formed by N-S-N triplets of saddles (S) and nodes (N) at 2 = 0. A snapshots 
of the skin friction field is shown. The solid lines between the nodes are thought as a guide to the eye. Saddles are always found 
in between two nodes close to this straight connection. The contours represent the temperature derivative dT/dz. 



FIG. 8: (Color online) Field lines of the surface vorticity field plotted together with the magnitude of the skin friction field, |s|, 
as background contour plot. The data are the same as in Fig. 
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FIG. 9: (Color online) Comparison of Lagrangian analysis (a-c) and corresponding Eulerian analysis (d-f). Three snapshots of 
the skin friction held dynamics are shown. The top row shows the attracting Lagrangian Coherent Structures as solid lines. 
The background shows the vertical temperature derivative. In the bottom row, held lines of the skin friction held and the 
node-saddle-node triplets are displayed. 


vorticity, uj±, for which uj± -3 = 0 holds. By comparing this figure and Fig. one can see that the unstable nodes 
of the skin friction field are surrounded by closed field lines of surface vorticity. The outermost closed field lines of 
the surface vorticity coincide with the local regions of enhanced shear which form between the regions where the cold 
(hot) plumes hit the bottom (top) plate and where new plume ridges are formed. This outermost closed field line is 
in some cases detected as a limiting streamline to which closed field lines of the surface vorticity converge from the in- 
and outside. This is observed for the upper left node in Fig. and is connected with field lines of the skin friction field 
that remain not perfectly radial further away from the critical point. 


D. Lagrangian analysis of line-like plume formation 

The Lagrangian analysis takes into account the time-dependent dynamics and is based on a set of 300 snapshots 
which are separated by dt = O.OlATf from each other. The data have been obtained at the same Ra, Pr and F. The 
integration time step is dt/20 and thus a total time of t = 4,.2Tf is spanned by this data set. The fine time resolution 
is necessary to conduct the Lagrangian analysis in order to determine the LCS. We integrate Lagrangian trajectories 
initialized on a fine grid by a fourth-order Runge-Kutta method and interpolate linearly with respect to time between 
the snapshots and by means of cubic splines with respect to the space. Approximations of the FTLEs are obtained by 
computing the maximum relative dispersion between neighboring trajectories, see e.g. [35]. The LCS are then extracted 
from the local maxima of the resulting scalar field. Figure 9 displays the result for three time instants ti. Here attracting 
LCS are obtained by a backward-in-time integration starting at tg.i = ti — Wdt. They are highlighted as the black 
solid lines in the top row contour plots of dT/dz. For comparison, we replot in the bottom row of the same figure the 
streamlines of the skin friction field together with the node-saddle-node triplets. The result confirms that the attracting 
LCS which are largely determined by the unstable manifolds of the saddles in the time-frozen fields are found at the 
line-like plume ridges. 
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VI. SUMMARY AND DISCUSSION 

We discussed the connection between the moving skeleton of critical points of the skin friction field and the vertical 
temperature derivative at the heated plate as a blueprint of the near-wall dynamics of line-like plume formation right 
above the heated plate. Node-saddle-node triplets are identihed as the typical configurations which form line-like plumes. 
Nodes are found where cold fluid hits the bottom plate. This is in contrast to former results in plane shear flows in 
which saddles and nodes form almost exclusively pairs rather than triplets. Beside the different geometry of the plates, 
there are several further differences between a boundary layer in a convection cell and a plane channel. For the RBC 
case, the unidirectional mean flow is absent. As it was discussed in Refs. [HUH!, the mean flow in a convection cell 
is fully three-dimensional and obeys its own complex slow dynamics (see also |39jf . Fluid turbulence is generated by 
buoyancy forces in the present case. Line-like thermal plumes, which form a skeleton that moves across the plate, are a 
typical feature that does not exist in a shear flow. The plume ridges enclose regions in which cold fluid is pushed back 
towards the bottom plate and generates unstable nodes of the skin friction held. Ring-like regions of high shear form 
around these unstable nodes. They can be identihed by closed held lines of the surface vorticity. 

We have also discussed that our observation does not violate the Poincare-Hopf theorem since x(M) = 1 for the 
present case. Finally, the application of Lagrangian Coherent Structure analysis showed that our ideas hold in a dynam¬ 
ical framework. The unstable manifolds of the saddles points visually coincide with attracting hyperbolic Lagrangian 
Coherent Structures. Attracting hyperbolic LCS are also found where the plume ridges are formed which subsequently 
rise into the bulk. 

With a view to the future work several points are interesting. First, the local how reversals which should be found in 
the vicinity of foci are rare as Fig. shows. This is similar to the channel how analysis. Very recently, they have also 
been found experimentally m in a shear how and are thought to play a central role for the formation of large-scale 
structures. Their role in RBC is yet open since the large-scale how is not uni-directional and depends on time. The 
percentage and dynamical relevance of foci could increase as the Rayleigh number grows. 

Second, we studied here one data set in order to demonstrate the ideas and concepts. It is certainly interesting 
and necessary to verify how robust the picture is for higher Rayleigh numbers. Furthermore, what happens when the 
Prandtl numbers are much larger or smaller than one and both boundary layers get decoupled? It can be expected that 
the dynamics is very diherent for the high-Prandtl case in comparison to the low-Prandtl case (at same Ra) simply due 
to the fact that the coherent large-scale motion is significantly weaker in the former case. Some of these analyses have 
been started very recently and will be reported elsewhere. 

Third, we saw that the skin friction fields at the top and bottom plates are not directly synchronized. The reason 
is that thermal plumes cannot rise straight to the top plate for the present Prandtl numbers. However, a large-scale 
circulation exists in a closed cell. For the present aspect ratio, a single large roll is expected [10] which couples the 
dynamics at the top and bottom plates. The particular way of coupling requires further detailed analysis. 



FIG. 10: (Color online) Direct numerical approximation of the regions around unstable nodes of the skin friction field. Dominant 
eigenvectors of a stochastic matrix sampled from Lagrangian trajectories highlight the dynamically distinct regions, which compare 
well with the corresponding results for the Lagrangian coherent structures in Fig. [^a). 


Finally, as especially the regions around unstable nodes are key to the plume formation, a direct and robust approx¬ 
imation is needed for further investigation. Set-oriented numerical methods using stochastic matrices are tailored for 


finding dynamically distinct regions in complex flows [miig. A preliminary study along this line is shown in Fig. HR 


The resulting partition of the bottom plate compares well with the corresponding LCS which is shown for the same data 
set in Fig. [^a). The application of such transfer-operator based approaches to the computational analysis of plume 
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formation in the full three-dimensional Rayleigh-Benard system is a subject for future research. In order to quantify for 
example the transport of heat and momentum which is associated with the sets that are shown in Fig. |10[ extensions 
of the existing mathematical framework will be necessary. 
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